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La presente invention se rapporte au domaine du 
traitement de 1' image. 

5 

La presente invention se rapporte plus 
particulidrement a un proc6d6 d' amelioration de 1' image 
d'empreintes digitales dans le but de faciliter 1' extraction 
des param£tres n6cessaires a une reconnaissance d'empreinte. 

10 

L'art ant^rieur connait d6ja par la demande 
WO03/079274 un proc6d6 d' amelioration des images 
d'empreintes digitales, oil les images enregistr6es sont 
trait6es a partir de blocs de 1' image initiale en utilisant 
15 des filtres de Gabor adaptes a des quantites 
caract£ristiques de chaque bloc. Plus precis&nent, le 
traitement decrit comprend les Stapes successives 
suivantes : 

- calcul des gradients locaux G de 1' image initiale / 
20 - calcul de l r espacement R entre les lignes, 

correspondant a une frequence ; 

- d6coupage de l 1 image de depart en diff brents blocs ; 

puis, pour chaque bloc : 
25 I 

- calcul de 1 1 orientation du bloc et de la variance de 

1 1 orientation a l'aide du gradient G ; j 

calcul de la frequence du bloc et sa variance & 
I'aide de R ; 

30 - construction de filtres de Gabor a partir de ces | 

param^tres ; 

- filtrage des blocs a I'aide des Gabor ; 

- ajout du bloc filtr£ a 1 9 image finale. 
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L' inconvenient d'une telle methode est que les 
orientations sont calcuiees a partir des gradients locaux de 
1' image initiale avant 1' application du ou des filtres de 
5 Gabor. Cette image pouvant etre de tr£s mauvaise qualite, 
les erreurs obtenues dans les calculs des gradients ont des 
consequences tr£s n6fastes sur les paramdtres choisis du ou 
des filtres de Gabor et done sur la quality de 1' image 
finale utilis6e pour 1' extraction. 

10 

I/art anterieur connait egalement par la publication 
de Hong et Al. « Fingerprint Enhancement » une methode 
d' amelioration d' empreintes . Cette publication mentionne 
l'int6r§t de filtrer 1' image par des filtres de Gabor avant 
15 de faire des choix d' orientations, Dependant, la publication 
utilise pour de tels filtres de Gabor une frequence moyenne 
fix6e ici & 60 cycles par hauteur d' image. 

L' inconvenient d'une telle methode est qu'en realite, 
20 il est fd&c de consid^rer que l'espacement des lign^S est 
fixe f et encore une fois, une telle hypothese est tr£s 
n6faste pour la qualite des r£sultats. 

On connait egalement, par la publication scientifique 
25 XP010614029 de Okombi-Diba et Al. « Segmentation of 
spatially variant image texture », un precede pour 
caract^riser les textures d' image et segmenter cette image 
selon ces textures. Le proc6d<§ se base sur 1' utilisation 
d'ondelettes de Gabor par des filtres de Gabor. Un tel 
30 proc6d£ est uniquement utilise pour la determination de 
textures . 
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La presente invention entend remedier aux 
inconvenients de l'art anterieur en calculant les 
orientations de 1' image apres un filtrage adapte a 1' image. 

5 Pour ce faire, la presente invention est du type 

d6crit ci-dessus et elie est remarquable, dans son acception 
la plus large, en ce qu'elle propose un proced6 
d' amelioration d'une image d'empreinte digitale comprenant 
au moins des etapes de : 
10 - normalisation de 1' image initiale ; 

determination de la zone utile de l r image / 

Et caract6rise en ce qu'il comprend en outre des 
6tapes successives de : 

decoupage de 1' image en une plurality de blocs ; 
calcul de la FFT de chaque bloc ; 

determination des frequences des blocs d' image & 
partir desdites FFT ; 

application de filtres de Gabor de param£tres 
determines par lesdites frequences ; 

determination des orientations de 1' image. 

Avantageusement, lesdits blocs comprennent des zones 
de chevauchement . 
25 

De preference, lesdites frequences des blocs sont 
determinees par la frequence la plus importante dans une 
bande de frequence de plus haute energie. 



15 



20 



30 De preference, ladite etape de determination des 

frequences des blocs d' image comprend en outre une etape 
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devaluation de la pertinence du calcul des valeurs desdites 
frequences . 

Avantageusement, en cas de non-pertinence d'une 
5 desdites frequences, ladite frequence est recalculee a 
partir d'une fonction de la FFT initiale. 

Selon un mode de realisation, la non-pertinence d'une 
desdites frequences est evalu£e par rapport a un seuil 
1 0 predetermine . 

De preference, ladite etape de determination des 
orientations de 1' image comprend en outre une etape 
devaluation de la pertinence du calcul des valeurs desdites 
15 orientations . 

Avantageusement, lesdites filtres de Gabor ont pour 
parametres d' orientation 0°, 22,5°, 45°, 67,5% 90°, 112, 5° f 
135°, et 157,5°. 

20 

De preference, ladite etape de determination des 
orientations comprend les etapes consistant a : 

reconstituer des images & partir desdits filtrages 
de Gabor desdits blocs ; 
25 - calculer l'intensite moyenne de chaque image 

filtree pour des zones de taille predeterminee ; 

cr6er une nouvelle image d' orientation contenant 
1' orientation du bloc de plus forte intensite ; 

cr6er une nouvelle image de qualite contenant 
30 l'intensite du bloc de plus forte intensite ; 

filtrer ladite image des orientations. 
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Selon un mode de realisation, le procede comprend en 
outre des etapes de creation d'une image finale fusionnee & 
partir desdites orientations, de binarisation et de 
squelettisation de ladite image finale fusionnee pour la 
5 detection des minuties. 

On comprendra mieux 1' invention & l'aide de la 
description, faite ci-apr&s a titre purement explicatif, 
d'un mode de realisation de 1' invention, en reference aux 
10 figures annex^es oCi : 

- la figure 1 illustre le pr6traitement de 1' image 
brute ; 

- la figure 2 illustre 1' amelioration de 1' image par 
15 filtrage. 

Illustre figure 1, le premier bloc fonctionnel selon 
1' invention est un bloc de pre-processing. II permet a 
1' image d'etre traitee ulterieurement dans des conditions 
20 satisfaisantes . " 

Ce bloc prend en entr6e 1' image brute lb issue du 
capteur d'empreintes digitales et lui applique diff6rents 
pretraiteraents. Ces pretraitements sont de type connu et 

25 consistent en une normalisation de 1' image ainsi qu f en une 
premiere determination de la zone utile de 1' image, c'est-a- 
dire la zone contenant 1' image d'empreinte. 

Pour la normalisation globale, on fixe la moyenne 
globale de tous les pixels de 1' image ainsi que leur ecart 

30 type cL une valeur pred6terminee, par exemple moyenne « 127 
et ecart-type=50 en niveau de gris. 
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Ainsi, si N=L*H est le nombre total de pixels de 
1' image brute et Ib(x) la valeur du xieme pixel de 1' image 
brute , on calcule : 

iV xd q 



Puis la deviation D b = ^^^Jb(xf-^ 



Et la nouvelle image normalisee est, si on a fix6 la 
nouvelle moyenne & Mo et la deviation k D 0 : 

10 

In(x)« M e i ^ v _ M * 

^6 



On applique 6galement un masque de type connu & partir 
des raoyermes et des 6cart-types locaux. A partir de cette 

15 zone utile, 1' image est retai!16e, Le retaillage de 1' image 
satisfait k deux objectifs : d'une part, elle permet de 
r6duire_J.a taille de 1' image afin de diminuer les temps de 
traitement lors de la g6n£ration du gabarit. Cette 
diminution de la taille de 1' image permet 6galement de 

20 limiter les risques de fausse detection de minutie a des 
endroits de 1' image ne faisant pas partie de 1'empreinte 
digitale elle-meme. D' autre part, elle permet de g6n6rer une 
image dont la taille sera un multiple de la taille d'un 
certain nombre de blocs se chevauchant, utilises pour 

25 1' amelioration de 1' image, afin de faciliter le traitement 
de ces blocs • 



30 



En sortie de cette premiere 6tape, on obtient une 
image normalisee retaill6e Ip. 
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A partir de cette image retaill^e Ip, on applique le 
bloc fonctionnel d' amelioration de 1' image lui-m§me illustr<§ 
figure 2. 

5 On decoupe d'abord 1' image en blocs de taille 

Tenh*Tenh, avec typiquement Tenh = 64 pixels. Ces blocs se 
chevauchent de Tchev afin d'6viter les effets de bord dans 
les traitements ult6rieurs. 

10 On applique ensuite un algorithme de FFT sur chacun 

des blocs afin d'obtenir leur transforra6e de Fourier. On 
peut noter que sur de tels blocs, 1' orientation des lignes 
ainsi que les espacements varient peu, done la valeur 
obtenue dans 1'espace reciproque est normalement tr£s 

15 localis£e. Cependant, il convient de d6finir une frequence 
moyenne pour le bloc. 

Cette frequence est obtenue de la fagon suivante : on 
applique d'abord l f algorithme permettant de calculer le 
spectre du bloc FFT en coordonn6es polaires. L' utilisation 

20 cfe coordonnees polaires est commode car elle permet de 
distinguer facilement les composantes du bloc representant 
la frequence de celles representant 1' orientation. On 
applique ensuite un algorithme de detection de la frequence 
la plus importante dans chaque bloc. Pour evaluer cette 

25 frequence, chaque bloc est decoupe en Nbandes bandes de 
frequence, avec typiquement Nbande «- 8. On choisit ensuite 
la frequence la plus importante dans la bande de frequence 
de plus haute §nergie. Cet algorithme permet aussi d'6valuer 
la pertinence de sa mesure. Cette valeur de pertinence est 

30 calcul^e en utilisant un rapport entre l'6nergie de la 
frequence la plus importante des bandes et l'6nergie moyenne 
des autres bandes. Pour certains blocs dont la mesure de 
frequence est jug6e non pertinente f une nouvelle Evaluation 
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est rdalisee en interpolant la valeur k partir de celle de 
ses voisins. 

Cet algorithme permet done d'obtenir une bonne 
5 estimation de la frequence moyenne des lignes d'empreintes 
pour chaque bloc, ainsi qu'un indice de confiance 
correspondant k la fid61it6 de cette estimation. 

De fagon conditionnelle, pour les blocs dont la 
10 quality d' estimation de la frequence est trop faible, un 
nouveau traitement est applique. La FFT est globalement 
multipli^e par elle-raeme afin d' augmenter les pics de 
frequence et d'eliminer les frequences faibles parasites. 
Une normalisation est ensuite effectuee afin de ne pas 
15 generer de donnees satur^es. Plus precisement, si ReFFT est 
la partie reelle de la FFT et ImFT sa partie 
imaginaire, on calcule pour chaque point, 

pseudoMagnitude{x) = (IteFFT(x) 2 + hnFFT(x) 2 ) n avec n-0 . 2 5 

typiquement et on choisit pour nouvelles coordonn^es de la 
*"20 FFT : ReFFT (x) = ReFFT (x) *pseudoMagnituae (x) et ImFFT(x) = 
ImFFT (x) * pseudoMagnitude (x) . La normalisation se fait 
alors en divisant chaque valeur complexe de la FFT par le 
Max des magnitudes du bloc. 

25 Le bloc fonctionnel d' amelioration est alors applique 

grace A cette carte des frequences Fr6q par des filtres de 
Gabor. 



Les filtres de Gabor sont des filtres selectifs en 
30 orientation et en frequence. lis permettent de faire 
ressortir des lignes espacees selon une certaine frequence 
et orient6es dans une certaine direction. Le principe de 
notre algorithme de filtrage est done d'utiliser un banc de 



WO 2005/071602 



9 



PCT/FR2005/050010 



f litres de Gabor, dont la frequence est d£termin£e par les 
frequences Fr6q definies plus haut. On utilise Nori 
orientations differentes pour definir les orientations du 
banc de filtre, avec typiqueraent Nori =8 pour des 
5 orientations de 0°, 22,5°, 45°, 67,5% 90% 112 f 5°, 135°, et 
157,5°. Les blocs sont filtr6s par les Nori f litres, et Nori 
images sont reconstitutes & partir des diff6rents filtrages 
des diff6rents blocs. 

10 A partir de ces Nori images, on procede £ la 

determination des orientations des images filtrees selon le 
processus suivant : on effectue une premiere estimation de 
1' orientation en calculant 1' intensity moyenne de chaque 
image filtr6e pour des blocs de taille TblocA, avec 

15 typiquement TblocA - 4* On ecrit alors dans une image Ori, 
1' orientation du bloc de plus forte intensite, et dans 
qual2, 1' intensity de ce bloc. 

La carte des orientations Ori ainsi obtenue pr6sentant 
20 de nombreuses discontinuites, >: on applique un filtre median 
afin de r6duire celles-ci et obtenir une representation plus 
homogdne . 

On recherche ensuite dans 1' image Ori ainsi filtree 
25 toutes les regions formees par des pixels connexes ayant la 
meme orientation. On va ensuite filtrer ces regions en 
fonction de leur taille et de leur orientation dans 1' image. 
Pour cela, on remplace les regions qui sont trop differentes 
de leurs voisins, selon un seuil pred£fini, par la valeur 
30 d' orientation la plus probable, calculee en fonction de 
1' orientation et de la quality des different s pixels voisins 
du pixel de la region trait6e. La carte de quality qua!2 est 
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mise & jour en fonction de l'intensite des pixels modifies 
lors de ce traitement. 



On applique ensuite le m§me traitement mais pixel par 
5 pixel, afin de r6duire les probl&nes dus a des changements 
brusques d' orientation. Pour effectuer ce traitement, on 
recherche tous les pixels dont la valeur est trop differente 
de celle de ses voisins. On applique alors apres une 
dilatation de 1 / image des pixels £ traiter le m£me 
10 algorithme de remplacement des orientations que pr6c6derament 
pour les pixels & traiter. La carte de qualite qual2 est 
mise d jour en fonction de l'intensite des pixels modifies 
lors de ce traitement 

15 on recherche ensuite une derniere fois les regions qui 

sont trop diff<§rentes de leurs voisines afin de determiner 
si les traitements effectu^s avant ont eu ou non un effet 
bdndfique. Si ce n'est pas le cas, on marque ces regions 
comme mauvaise dans la carte de quality qual2. 

20 

Enfin, on applique un filtre median 4 cette carte de 
quality - 

En fonction de la carte des orientations calculees 
25 pr6cedemment, on s61ectionne les pixels a ecrire dans 
1' image fusionnee finale. Ainsi, en fonction de la carte des 
orientations, on choisit chaque pixel dans une des Nori 
images filtrdes. Afin de rendre les transitions plus douces, 
une interpolation entre les diff§rentes orientations est 
30 rdalisee aux fronti&res de chaque bloc. 



L' image fusionnee ainsi obtenue est alors binarisee et 
squelettisee de fagon connue pour permettre la detection des 
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mimities f la carte de quality qua!2 servant & determiner la 
pertinence de chaque minutie. 

L' invention est decrite dans ce qui precede & titre 
d'exemple. II est entendu que l'homme du metier est a m§me 
de r£aliser diff^reiites variantes de l f invention sans pour 
autant sortir du cadre du brevet. 
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REVENDICATIONS 

1. Proc6de d' amelioration d'une image d'empreinte 
digitale comprenant au moins des etapes de : 

5 - normalisation de 1' image initiale ; 

determination de la zone utile de 1' image ; 
et caracterise en ce qu'il comprend en outre des 
etapes successives de : 

decoupage de 1' image en une pluralite de blocs ; 
10 - calcul de la FFT de chaque bloc ; 

determination des frequences des blocs d' image/ <S. 
partir de ladite FFT ; 

application de filtres de Gabor de paramdtres 
determines par lesdites frequences ; 
15 - determination des orientations de 1' image A partijc 

des images filtr6es issues de 1' etape precedente. 

2. Procede d f amelioration d'une image d'empreinte 
digitale selon la revendication 1, caracteris6 en ce que 

20 lesdits blocs comprennent des zones de chevauchement . 

3. Procede d f amelioration d'une image d'empreinte 
digitale selon la revendication 1, caracterise en ce que 
lesdites frequences des blocs sont determinees par la. 

25 frequence la plus importante dans une bande de frequence de 
plus haute energie. 

4. Procede d' amelioration d'une image d'empreinte 
digitale selon la revendication 1, caracterise en ce que 

30 ladite etape de determination des frequences des blocs 
d' image comprend en outre une etape devaluation de la 
pertinence du calcul des valeurs desdites frequences. 
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5. Procede d' amelioration d'une image d'empreinte 
digitale selon la revendication 1, caract6rise en ce que 
ladite etape de determination des orientations de JL' image 
comprend en outre une etape d' evaluation de la pertinence du 
calcul des valeurs desdites orientations. 

6. Proc6de d r amelioration d'une image d'erapreinte 
digitale selon la revendication 4, caracterise en ce qu'en 
cas de non-pertinence d f une desdites frequences, ladite 
frequence est recalcuiee & partir d'une fonction de la FFT 
initiale. 

7 . Proced6 d' amelioration d' une image d' erapreinte 
digitale selon la revendication 6, caract6rise en ce que la 
non-pertinence d'une desdites frequences est evalu«§e par 
rapport & un seuil predetermine. 

8. Procede d' amelioration d'une image d'empreinte 
digitale selon la revendication 1, caracterise en ce que 
lesdites filtres de Gabor ont pour param^tres d' orientation 
0% 22,5% 45% 67,5% 90% 112,5% 135% et 157,5% 

9. Procede d' amelioration d'une image d'empreinte 
digitale selon la revendication 1, caracterise en ce que 
ladite etape de determination des orientations comprend les 
etapes consistant k : 

reconstituer des images a partir desdits filtrages 
de Gabor desdits blocs ; 

calculer l'intensite moyenne de chaque image 
filtree pour des zones de taille predeterminee ; 

creer une nouvelle image d' orientation contenant 
1' orientation du bloc de plus forte intensite ; 
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cr6er une nouvelle image de qualite contenant 
1' intensity du bloc de plus forte intensite ; 

filtrer ladite image des orientations. 

10. Proc6d6 d' amelioration d'une image d'empreinte 
digitale selon la revendication 1, caract&ris6 en ce qu'il 
comprend en outre des etapes de creation d'une image finale 
fusionnee i partir desdites orientations, de binarisation 
et de squelettisation de ladite image finale fusionnee. 
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